A new universality class of aging in glassy systems 
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Glassy systems are ubiquitous in nature, and are characterized by slow relaxations to equilibrium 
without a typical timescale, aging and memory effects. Understanding these is a long-standing 
problem in physics. We study the aging of the electron glass, a system showing remarkable slow 
relaxations in the conductance. We find that the broad distribution of relaxation rates leads to a 
universal relaxation of the form log(l + t w /t) for the usual aging protocol, where t w is the length of 
time the perturbation driving the system out-of-equilibrium was on, and t the time of measurement. 
These results agree well with several experiments performed on different glasses, and examining 
different physical observables, for times ranging from seconds to several hours. The suggested 
theoretical framework appears to offer a paradigm for aging in a broad class of glassy materials. 



PACS numbers: 71.23.Cq, 73.50.-h, 72.20.Ec 

Aging is one of the most distinct characteristics of 
glasses. Crudely speaking, it is the dependence of the 
form of the relaxation to equilibrium on the time an 
external perturbation was acting on the system. Aging 
and slow relaxations have been experimentally observed 
in various systems. These range from spin-glass to 
structural glasses @, Q-gels vortices in superconduc- 
tors ferroelectrics [y], colloids and granular systems 
0| and the electron-glass system Analytic models 
and simulations showing aging have been of considerate 
interest in the past decades [9|. In general, the return to 
equilibrium can be characterized by a function F(t,t w ), 
where t w , the 'waiting-time', is the time the perturbation 
was on (during which the system is being 'aged'), and t 
is the time from switching off the perturbation until the 
measurement. We refer to the complete protocol as an 
'aging experiment'. Full (or simple) aging is said to hold 
when F{t, t w ) = f{t w /t). 

Experiments performed on the electron glass system, 
measuring the conductance as a function of time when 
the system is perturbed, show remarkable data collapse 
when time is measured in units of t w Q ■ This proves ex- 
perimentally that full aging is obeyed. For times t <C t w 
a logarithmic relaxation is observed. At times t of or- 
der t w deviations from the logarithm are seen. Sim- 
ilar results were obtained for granular aluminum (loj 
and structural glasses 0] which suggest that the phys- 
ical mechanisms leading to such logarithmic relaxations 
may be ubiquitous. In this paper we suggest a model 
which should be applicable for several types of physical 
realizations, and yields the result that full aging indeed 
holds. Moreover, it gives under certain simplifying as- 
sumptions f(x) ~ log(l + a;), which fits the experimental 
data well for several decades in time. 

The Model- To start with, let us consider N localized 
disordered states and M < N interacting electrons, with 
a coupling between the electrons and a phonon reser- 
voir llj. After ensemble averaging over the phonon 



reservoir, the averaged occupation numbers can take any 
value between and 1, and their dynamics is described 
by a set of rate equations: ^£ = Sj^i ~ 7i,j'i where 
due to the Coulomb interaction and the Pauli principle 
the rates ji j, which describe the phonon-assisted transi- 
tions between electronic states, depend on the set of oc- 
cupation numbers {n-i}. Linearization close to (one of the 
many) metastable states leads to an equation = A5n, 
where Sn is a vector of the deviations from the metastable 
state, and A a matrix depending on the properties of 
the metastable state. Ref. ll| explains in detail the 



procedure, and shows that the distribution of eigenval- 
ues of the matrix is given approximately by P(A) ~ 1/A 
within a large frequency window, essentially arising from 
the exponential dependence of the matrix elements on 
the physical parameters (distance and energy difference). 
The 1/A distribution leads to a logarithmic decay in time, 
as can be seen by Laplace transforming it. 

From now on we consider quite generally a system 
whose dynamics is described by a linear relation = 
A ■ 6n, and with the eigenvalue distribution of the matrix 
A described approximately by P(A) ~ 1/A. The mean- 
field approximation leads to this result starting from a 
microscopic Hamiltonian, but we believe that the results 
should apply for a broader range of models. At the end 
of the paper we demonstrate this by applying the theory 
to the experimental data of Ref. [3(] . 

We will assume that the system contains another exter- 
nal control parameter we shall call V g , such as the gate 
voltage in the electron glass experiment Q. Changing 
an external system parameter will shift the metastable 
states' positions. Thus, upon small changes in V g , the 
linearized matrix A will slightly change. Let us consider 
the system in two metastable states a and b correspond- 
ing to values V a and Vj, respectively. We shall denote 
the linearized matrix A in the state a(b) by A a {A b ), 
and the eigenvectors a q (b q ). The metastable state a is 
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characterized by a configuration ft a , and the state & by a 
configuration n& = ri a + An. 

The experimental protocol is as follows: 

I The gate voltage is V a . The system is assumed to 
'equilibrate' to a mctastable state. 

II The gate voltage is changed to Vb = V a + SV, the 
system 'ages' for a time t w : the system approaches the 
new equilibrium b. 

III The gate voltage is changed back to a value V a . 
The system is in the process of returning to the original 
equilibrium a. 

The conductance is measured throughout the exper- 
iment. The measurement time is much shorter than 
the relevant relaxation times. For convenience we de- 
fine t = at the end of stage II, when the gate voltage is 
changed to its initial value V a . 

During stage I, the system reaches, presumably, one of 
the many metastable states (a), and lingers there. At the 
beginning of stage II (t = —t w ) we change the potential 
landscape by shifting suddenly the gate voltage from V a 
to Vb- The occupations at this time are still equal to 
the mctastable state occupations ri a . During stage II it 
starts, however, to relax to the new metastable state b. 
The new metastable state b is not far (in phase space) 
from the original one a , but since the relaxation to it 
is composed of a broad range of timescales, the system 
takes a long time to fully equilibrate to state b. At the 
beginning of stage III (t = 0) the gate voltage is changed 
back to its initial value. During stage III, the system is 
relaxing back to a, the metastable state it was initially 
in [12J ■ This process is described schematically in Fig. [IJ 
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along with our theoretical predictions. 

Let us proceed to the calculation. In stage I, the sys- 
tem is in the metastable configuration n a . At the begin- 
ning of stage II, the system is in the state n a = fib — An. 
Working in the eigenbasis of b, —An = ^2 q c q b q . This 
equation determines the decomposition weights c q (which 
are not simply the projections onto the eigenvectors, 
since the matrix A b is not hermitian [HI). Therefore 
during stage II (the aging process), the configuration of 
the system is described by: 




ft(t) = n b + 



J2 c i^i e ~ 



x «( t+t ™\-t w <t<0, (1) 



FIG. 1: Schematic demonstration of the aging process. Ini- 
tially, the system is in a metastable configuration a. When 
the landscape changes due to a change of gate voltage, the 
system relaxes towards the new minimum b, during time t w - 
When the landscape returns to its initial form, the system 
relaxes to a. Below we schematically show our results, for the 
time dependence of the distance to the current metastable 
state (a or 6) Sn(t). At first, the system sits at a metastable 
state. At a certain time (-200 sec. for the figure) the system 
parameter (i.e., gate voltage) is changed, and Sn(t) jumps and 
then relaxes logarithmically. At time t w the system parame- 
ter is changed to its initial value, from which the relaxation 
takes the form log(l + t w /i) (Eq. §2$). The inset shows the 
experimental result for a typical electron glass aging experi- 
ment (courtesy of Z. Ovadyahu). The lower cutoff used in the 
theory, relevant only for stage II, is a fitting parameter, and 
was taken to be 50 hours. 



where A^ are the relaxation rates (the eigenvalues of 
the linearized matrices A a and A b are real and nega- 
tive, corresponding to pure decay). Since the eigen- 



vectors are localized ll|, we can assume that c q are 



random variables whose distributions and typical values 
do not depend explicitly on q. Thus, for a large sys- 



tern: | £ 



q c q b q e 



-dX = Ei(\ max At)- 



Ei(XminAt), where Ei(x) is the exponential integral 
function. For times At — ► the integral approaches a 
constant, which by construction is |An|, while for times 



At ^> 1/A m j„ it falls off exponentially. The interest- 
ing behavior occurs for the large time window l/X ma x < 
At < l/X m in, where we can approximate the integral by 
7b— log[AtA m i n ], with is the Euler constant. This im- 
plies that during stage II we should observe a logarithmic 
relaxation, which is indeed observed experimentally [l3| . 

In stage III, it is convenient to work in the eigenvec- 
tor basis of metastable state a. We can express the 
configuration of the system at the end of stage II as 
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ft(t = 0) = n a + Ari+J2 q Cqb q e x i tw . Since the perturba- 
tion is small, to lowest order we can replace b q by a q and 
substitute the eigenvector decomposition of An, obtain- 
ing: ft(0) = ri a + J2 q c qa q {e" x i tx " — 1). Note that during 
the waiting-time, each of the eigenmodes approaches b 
with a characteristic time 1/A g . Therefore only modes 
obeying Xt w > 1 had enough time to relax to b. When 
the voltage is changed to its initial value, this relaxation 
to b is in fact an excitation with respect to the metastablc 
state a. Thus, slowly relaxing modes are also slow to be 
excited, which is the crux of the matter. The relaxation 
back to equilibrium during the last stage of the protocol 



-l)e 



is described as: n(t) = ri a + ^2 q c q a q ( 

We can approximate A° ~ A 
occupancy vector at time t: Sn(t) = ^2„c q a q ( 



A g , and obtain for the 



-x t u 



l)e qt . Defining 6n(t) = \Sn(t)\, we obtain: 



Sn(t) = C[~fE-log(t+t w )}-C[~/ E -log(t)] 



where C = 



[Agj 



C\og(l+t w /t), 
(2) 

is a non- universal constant, 



log(A max /A 

which depends on the system parameters. 

Notice that is a function of t/t w , correspond- 

ing to full aging. This is due to the fact that for 
f(x) ~ log(x), the difference f(t + t w ) — f(t) is only 
a function of the ratio t/t w . It can be proven that only 
a logarithmic function has this characteristic. Since f{x) 
was essentially a Laplace transform of the eigenvalue dis- 
tribution, this shows a remarkable relation between the 
1/A distribution and full aging. Fig. Q] summarizes the 
theoretically expected time dependence throughout the 
experiment, and shows the results of a typical electron 
glass aging experiment. 

So far we have discussed the aging of the average occu- 
pation numbers. To relate this to the measured observ- 
able in the electron glass experiments, the conductance, 
we need to know the relation 5o~{5n). We emphasize that 
for any such relation we would obtain full aging, since Sn 
has full aging. In the simplest approximation, one would 
assume a linear relation between the two. In other words 
Sn(t), essentially the distance in phase space to the lo- 
cal minimum (the metastable state), is proportional to 
the excess conductivity. This may arise due to various 
mechanisms. One may intuitively regard perturbation in 
occupation number as an effective excess temperature, 
raising the conductivity. A more subtle mechanism would 
arise due to slowly relaxing modes, influencing the con- 
ductivity of the system via the long-range Coulomb in- 
teractions [ill. [T3|. This is an important point, that may 
distinguish the electron glass system from spin glasses, 
and we intend to elaborate on it in future work. 

Comparison to experiment.- Fig. [2] compares the the- 
oretical form in stage III and the experimental data of 
Ovadyahu's group [8(. The experimental data shows a 
logarithmic regime for times shorter than t w , and power- 



law regime, for times much longer than t w . These two 
limits, as well as the crossover between them, are con- 
tained in the theoretical expression. The significance 
of the fit shown in the figure is the linear dependence 
of the excess conductivity on the predicted function 
log(l + t w /t). There are no fitting parameters other than 
the proportionality constant, which is non-universal. 

The logarithmic behavior predicted for short times is 
an important distinction between the theory we present 
and various theories predicting a stretched exponential 

behavior, e^T & At short times, compared to t w , 
the stretched exponential reduces to a power-law, which 
is inconsistent with the experimental data, as was demon- 
strated in (l3j |. 

For even longer times (on the scale set by t w ), devia- 
tions from the theory occur. These can be corrected by 
taking a non- linear relation between Sn and 5a. In order 
to fit correctly the experimental data, the corrected re- 
lation must be sublincar. An additional sublincar term 
proportional to y/Sn provides excellent correspondence to 
the experimental data (not shown). That such a sublin- 
ear relation may exist is supported by a different exper- 
iment, where the conductance is measured as a function 
of the gate voltage in a relatively fast scan 

Generality of the theory.- It is possible that this model 
is in fact only a caricature of the real scenario: possibly, 
the system explores many metastable states (local min- 
ima of the potential landscape). If we consider the vector 
of the probabilities to remain in each state, the transi- 
tion between different minima are also exponentials of 
a broadly distributed quantity [Til . This implies that 
the 'ingredients' used in our calculation are still valid: 
if we consider the vector of probabilities of the system 
to be in each of the minima, the (classical) dynamics of 
the system will be described by a master equation of the 
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FIG. 2: The excess conductance is plotted as a function 
of the theoretical expression log(l + t w /t) . There are no 
fitting parameters in this expression. The data shown is for 5 
experiments, corresponding to different colors, with different 
waiting times t w of 20, 60, 180, 540 and 1620 seconds. The 
ratio t/tw varies from 0.0025 to 36. We thank Z. Ovadyahu 
for providing us with the data. 
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type ^£ = AP, where P is a vector of probabilities to 
be at different minima. A wide, slowly varying, distri- 
bution of barriers would again lead to a 1/A distribution 
of eigenvalues of the matrix A, leading to the relaxation 
described by Eq. ([2]). 

We demonstrate the generality of the theory by an- 
alyzing a different experiment. Ludwig et al. [3] mea- 
sured the dielectric constant of Mylar when pushed out- 
of-equilibrium, using a similar aging protocol to the one 
described. We consider the data from the first mea- 
surement in each of their aging experiments, for which 
the system is fully equilibrated. The data presented in 
the first reference of [3[ shows that for times short com- 
pared to t w , the signal is proportional to \og(ct w /t), with 
c = 1.02 extracted from the data. This fits well with our 
prediction of log(l + t w /t). We emphasize that in con- 
trast to our theory, the theory used in Ref. [3] gives a 
non-universal value of c, depending on the system pa- 
rameters. 

Summary.- We have shown how full aging may nat- 
urally occur in a model for the electron glass, starting 
with a microscopic picture, and using a master equa- 
tion to analyze the dynamics of the system. Under sim- 
plifying assumptions, we have shown that the deviation 
from the metastable state should follow a universal form 
~ \og(l + t w /t), where t w is the length of time the pertur- 
bation driving the system out-of-equilibrium was on, and 
t the time of measurement. This agrees well with experi- 
mental data of the conductance relaxation in the electron 
glass InO, within a large time window, ranging from sec- 
onds to several hours. The theoretical predictions also ex- 
plain experiments performed on a totally different glass, 
the plastic Mylar, where the dielectric constant is mea- 
sured. Although the system is completely different, the 
universality of the underlying distribution of relaxation 
rates P{X) ~ 1/A lead to our general result for the relax- 
ation. This suggests that our model of relaxation arising 
from a broad 1/A spectrum may be a paradigm for aging 
in various glasses. 

We thank Z. Ovadyahu and S. Ludwig for kindly pro- 
viding us with their experimental data, as well as for 
illuminating discussions. 
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